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High-power and highly directional semiconductor cylinder-lasers based on an optical resonator with 
deformed cross section are reported. In the favorable directions of the far-field, a power increase of up 
to three orders of magnitude over the conventional circularly symmetric lasers was obtained. A "bow- 
tie"-shaped resonance is responsible for the improved performance of the lasers in the higher range of 
deformations, in contrast to "whispering-gallery" -type modes of circular and weakly deformed lasers. 
This resonator design, although demonstrated here in midinfrared quantum-cascade lasers, should be 
applicable to any laser based on semiconductors or other high-refractive index materials. 
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Lasers consist of two basic components. First, the ac- 
tive material in which light of a certain wavelength range 
is generated from an external energy source, such as elec- 
tric current; second, the laser resonator, which contains 
the active material, provides feedback for the stimulated 
emission of light. The resonator largely influences the 
special features of the emitted light: power, beam direc- 
tionality, and spectral properties, as well as the laser's 
physical features such as size and shape. Semiconduc- 
tor lasers are the most widely used and versatile class 
of lasers. Their most common resonators are Fabry- 
Perot cavities, in which two cleaved semiconductor crys- 
tal planes act as parallel mirrors, reflecting the light back 
and forth through the active material. 

There have been many attempts to improve resonator 
properties. In particular, an increase of the reflectivity 
of the resonator mirrors is highly desirable. This allows 
low thresholds for the onset of laser action and a smaller 
volume of active material with concomitant moderate en- 
ergy requirements and the ability to pack the lasers in a 
small space. 
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One excellent example is the development of micro- 
disk semiconductor lasers (1). These lasers exploit total 
internal reflection of light to achieve a mirror reflectivity 
near unity. Micro-disk, -cylinder or -droplet lasers form a 
class of lasers based on circularly symmetric resonators, 
which lase on "whispering-gallery modes" of the electro- 
magnetic field (2,3,4). In such a mode light circulates 
around the curved inner boundary of the resonator, re- 
flecting from the walls of the resonator with an angle of 
incidence always greater than the critical angle for total 
internal reflection, thus remaining trapped inside the res- 
onator. There are only minute losses of light caused by 
evanescent leakage (tunneling) and scattering from sur- 
face roughness. This principle allowed the fabrication 
of the world's smallest lasers (2). Besides possible fu- 
ture applications in optical computing and networking, 
micro-lasers are of strong interest for research problems 
of cavity quantum electrodynamics, such as resonator- 
enhanced spontaneous emission and threshold-less lasers 
(5). Small resonators may also serve as model systems 
for the study of wave phenomena in mesoscopic systems, 
particularly in the regime where motion is fully or par- 
tially chaotic. Recent examples are the quantum me- 
chanics of electrons confined in asymmetric "boxes" , such 
as quantum-dots, stadia, and quantum corrals (6), and 
asymmetric microwave cavities with their strong connec- 
tion to quantum chaos theory (7). 

However, as a serious disadvantage, the tiny 
"whispering-gallery" -type lasers lack high output power 
and directional emission because of the high reflectivity 
mirrors and the circular symmetry. Attempts to improve 
this deficiency by making gratings or small indentations 
on the circumference are so far not very promising (8,9). 
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Figure 1: Scanning electron microscope image of the side and 
top view of a flattened quadrupolar shaped cylinder laser. The 
quadrupolar deformation parameter is f K 0.16. Side view: The 
laser waveguide and active material is entirely contained in the disk 
with vertical side walls, total thickness d = 5.39 fim sitting on a 
sloped Indium-Phosphide pedestal. Light emission occurs in the 
plane of the disk. Top view: The top face of the laser is shown in 
medium gray, the electric contact in light gray. The laser bound- 
ary is very well described by an exact flattened quadrupole with 
e = 0.16 which is drawn as dashed red line over the circumference. 
The shape of the top electrical contact approximately parallels the 
laser shape. 

We now show in experiment and theory how a new res- 
onator design which incorporates chaotic ray motion can 
dramatically increase the output power and directionality 
of such lasers. This effect is demonstrated in semiconduc- 
tor quantum-cascade lasers emitting in the mid-infrared 
wavelength region (10). 

Recent theoretical work has provided insight into 
the behavior of "asymmetric resonant cavities" (ARCs), 
"whispering-gallery" resonators with smooth deforma- 
tions from cylindrical or spherical symmetry (11-14). 
The ray dynamics in these deformed resonators is known 
to be either partially or fully chaotic in the generic 
case (13). The most well-studied example is a two- 
dimensional (2D) resonator with a quadrupolar deforma- 
tion of the circular boundary, described in polar coordi- 
nates (r, (f>) by r(<f>) oc (1 + e cos(2 </>)), where e is the 
deformation parameter. Partially chaotic "whispering- 
gallery" modes in these resonators have shown directional 
lasing emission in low-index materials (index of refrac- 
tion n < 2, such as glass fibers or cylindrical dye jets) 
(12). The origin of the directional emission is the fol- 
lowing (11): the deformed boundary causes the angle of 
incidence of a ray in a "whispering-gallery" mode to fluc- 
tuate in time. Eventually a ray trapped by total internal 
reflection impinges on the boundary below the critical 
angle and escapes by refraction. However, it was not 
appreciated in this earlier theoretical work that, in high- 
index materials, qualitatively different modes not of the 
"whispering gallery" type might be relevant to the lasing 



properties. 

Here we focus on semiconductor lasers that have an ef- 
fective index of n w 3.3 and a deformation of the bound- 
ary best described by r(</>) oc y/\ + 2 e cos(2 <p), which we 
will refer to as a "flattened" quadrupole. In general, one 
can parameterize the boundary of any convex resonator 
by an arbitrary Fourier series for r(cj>); the above param- 
eterization is chosen for convenience, because it is simply 
analyzed and describes the actual resonator shapes quite 
well (see Fig. 1). We show that for small deformations 
e, the basic picture of chaotic "whispering gallery" or- 
bits escaping refractively, as described above, still holds 
for the high-index semiconductor material. However, we 
also present strong experimental evidence that at larger 
deformations a different type of laser resonance emerges 
and is responsible for highly directional and high power 
emission. Unlike the chaotic "whispering gallery" modes 
of smaller deformations, these so-called "bow-tie" reso- 
nances are stable resonator modes surrounded on all sides 
(in phase space) by chaotic motion. 

This new class of laser-resonators, based on the smooth 
deformation of a regularly shaped monolithic cavity, is 
universally applicable to semiconductor lasers or solid- 
state lasers based on high refractive index material. It 
is to some extent related to ring-lasers with resonators 
formed by assembly of several distinct flat or curved mir- 
ror surfaces, such as the Ti:Sapphire laser (15) or opti- 
cally pumped monolithic solid-state lasers (16). Never- 
theless we chose the quantum cascade (QC) laser, as it 
is particularly suited for 2D ' 'whispering-gallery" geome- 
tries as shown by a recent work on QC microdisk laser 
(17). It is based on a transition between quantized con- 
duction band states of a cascaded InGaAs/InAlAs cou- 
pled quantum- well structure (intersubband transition). 
As such, the selection rule of the optical transition al- 
lows light emission only in the 2D plane with polariza- 
tion normal to the quantum well layers; that is transverse 
magnetic (TM) polarization (18). Therefore, no light is 
lost vertical to the laser plane. Furthermore, the QC 
laser is a unipolar device based only on electron trans- 
port, unlike diode lasers. Thus, in contrast to most con- 
ventional semiconductor lasers, the surface can not cause 
excess unwanted non-radiative recombination of electrons 
and holes. Finally, the wavelength is comparatively large 
(several micrometers) and the material used is the well- 
understood InGaAs/InAlAs system. This choice reduces 
the importance of roughness (Rayleigh) scattering and 
makes it easier to fabricate complex shapes. 

Device structure and experimental procedure. 
The lasers are slab-waveguide structures made from a 
Gao.47Ino.53As/Alo.4gIno.52As heterostructure grown by 
molecular beam epitaxy (MBE) on InP substrate. The 
waveguide core contains the QC laser active material, 
designed to emit light of A = 5.2 /xm wavelength. This 
active material has been used previously (19) and can 
be considered a mature and optimized design for high 
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quality laser performance. 

The waveguide core is sandwiched between two 
cladding layers (20-22). The entire waveguide is designed 
to be symmetric and such that the lasing mode (the low- 
est order TM mode) has almost no (< 0.5%) overlap with 
the InP substrate. This design prevents possible detri- 
mental effects from light coupling into the substrate. 

The cylinder lasers are fabricated by optical lithog- 
raphy and wet chemical etching. The quadrupolar-likc 
shape is obtained starting from a resist pattern that is 
composed of two semi-circles connected by a rectangle 
(stadium-shape). The samples are then etched until deep 
mesas are obtained (23). Because of the smoothing ac- 
tion of the etchant the straight section of the etch mask 
bends toward the curved parts rendering a quadrupole- 
like shape. Figure 1 shows the top and side view of a 
laser with deformation e w 0.16. 

The top view shows that the edge of the resonator fol- 
lows very well the shape of an exact flattened quadrupole. 
Ohmic contacts are applied to the front and back sur- 
face of the lasers by using non- alloyed Ti/Au and 
Ge/Au/Ag/Au, respectively. 

Several sets of samples were fabricated. The deforma- 
tion parameter e was varied in 10 steps from to « 0.2. 
Two different sizes were investigated in order to quantify 
and rule out size-dependent effects, one with the short 
diameter w 50 /im. and the long diameter varying from 
s=s 50 jum (e = 0) to ~ 80 /im (e w 0.2), the second with 
the short diameter w 30 ^m and the long diameter vary- 
ing from « 30 /an to w 50 /im (24). The measurements 
described below showed that effects arising from the in- 
crease in cavity cross-section with increasing e (by less 
than a factor of 2, for < e < 0.2) are negligible com- 
pared with those introduced by the deformation. 

The measurements were mainly performed by contact- 
ing the individual cylinder laser with a micro-probe in 
a cryogenically cooled micro-positioner stage. To obtain 
the far-field pattern, the individual laser was mounted on 
a sample-holder that was rotated inside the probe stage. 
The lasers were driven with current pulses (duration 50 
ns, repetition rate « 40 kHz), and the light-output was 
measured with a cooled HgCdTe-detector and a lock-in 
technique. To improve power output and avoid excess 
current heating the data presented here were taken at 40 
to 100 K heat sink temperature. Nevertheless the max- 
imum pulsed operating temperature of the lasers is 270 
K. The spectral properties were measured using a Fourier 
transform infrared (FTIR) spectrometer. 

The lasers emit light according to their symmetry into 
all quadrants of the 2D laser plane. The experimental 
set-up did not allow to measure the spatially integrated 
power. Therefore we collected the laser output into an 
appropriate aperture. Its center angle is varied for the ac- 
quisition of the far-field pattern. The light output verti- 
cal to the laser plane is broadened by diffraction and was 
measured integrated over the vertical extension. We in- 



troduce a polar coordinate system (r, <f>) such that = 0° 
indicates the direction along the elongated (major) axis. 
Accordingly, <f) — 90° denotes the direction of the com- 
pressed (minor) axis. Hence a measurement taken at 
4> = 0° has the detector facing one point of highest cur- 
vature of the deformed laser. 

Power output and beam directionality. The de- 
formed cylinder lasers provided both a dramatic increase 
of the emitted power and directionality (Figs. 2 and 3). 

Light output measurements for various lasers as a func- 
tion of their deformation e are shown in Fig. 2A. To gen- 
erate this plot, the maximum obtainable peak power was 
recorded for each laser by optimizing the pulsed drive 
current (25). The collecting aperture (slit width cor- 
responding to 15°) was oriented around <j> = 0°. Note 
that this set-up precludes the observation of any changes 
in the far-field directionality with deformation. Similar 
measurements were performed with the aperture oriented 
around <j> = 45° and <j> = 90°. The striking result is the 
strong (quasi-exponential) increase of the collected opti- 
cal power with deformation. For the largest deformation 
under consideration (e w 0.2), a power increase of a fac- 
tor of w 50 with respect to the circular case is observed. 
Figure 2A shows a representative measurement taken at 
= 0°. 

The absolute output power was measured in some 
highly deformed laser devices by bonding and mounting 
them in a calibrated set-up usually used with Fabry-Perot 
type lasers. One example, obtained from a laser with 
e « 0.2, is shown in Fig. 2B. The light-collecting aper- 
ture was increased to its maximum size, and the sample 
was tilted to detect roughly the optical power in an an- 
gle from +40° to +100°. The choice of this aperture, 
which exploits far-field anisotropy, will become clear be- 
low. A peak output power of ~ 10 mW at 100 K was ob- 
tained. This value is approximately three orders of mag- 
nitude greater than that obtained from the non-deformed 
(circular cylindrical) laser or previous conventional cir- 
cular QC-disk lasers (17). For a weakly deformed laser 
(e = 0.06) we estimate a peak power output of « 50 //W 
(when measured with comparable collection efficiency as 
the laser of Fig. 2B). 

A quasi-exponential increase of the collected power 
with increasing deformation (similar to the one shown 
in Fig. 2A) has been measured in numerous sets of lasers 
of various - flattened and less flattened - quadrupolar 
shapes and sizes, and with various orientations of the 
aperture. Thus it appears that the power increase is a 
reliable, universal effect. However, the increase in out- 
put power per unit angle is rather closely entangled with 
the actual variation of the far-field pattern with defor- 
mation. In fact, in our lasers the power increase with 
deformation results from the lasing of different types of 
modes in different ranges of e. There is a cross-over at 
intermediate deformations (e « 0.12) from emission via 
"whispering-gallery" modes, which dominates at smaller 
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Figure 2: (A) Maximum peak optical power from various lasers as 
function of their quadrupolar deformation parameter e. The aper- 
ture with width 15° is centered around zero degree. Data from 
two independent sets of lasers are given. The power output is nor- 
malized to the power of the respective circular cylinder laser. The 
scatter of the data is due to the coarseness of the mode spectra 
(single or few versus multiple modes lasing). (B) Light output- 
to-current characteristics of a quadrupolar cylinder laser with de- 
formation e tss 0.2. The collecting aperture ranges from +40° to 
+100° (the polar coordinate system is described in the main text 
and in Fig. 3C). The kink around 400 mA indicates the onset of a 
second lasing mode. The measurement has been performed at 100 
K heat sink temperature. The lasers have been tested up to 270 K. 

deformations, to laser emission from "bow-tie" modes 
which do not exist below e rs 0.10 but dominate the high 
deformation regime. 

In addition to the strong increase in power output, 
the deformed lasers can also provide strong directional- 
ity. The results of the far-field measurements are summa- 
rized in Fig. 3A and 3C. As expected, the circular cylin- 
der laser displays no directionality of the emission. At 
small deformations (e < 0.10), the far- field is only weakly 
structured with an increased emission in direction of the 
minor axis compared to the major axis. Figure 3A shows 
the increase of the output power with e, collected around 
0° and 90°. Both curves rise approximately exponen- 



tially, as discussed in the previous section, but "faster" 
for 4> — 90°; in this case the exponent is increased by a 
factor of w 2 with respect to the <f> — 0° case (26). 

This observation is consistent with the expected be- 
havior of deformed "whispering gallery" modes with an 
average angle of incidence near the critical angle defined 
by sin(x c ) = 1/n, where n — 3.3 is the effective refrac- 
tive index of the laser waveguide. At zero deformation 
such a mode has a conserved angle of incidence and emits 
isotropically and uniformly via evanescent leakage from 
all points at the boundary (neglecting disorder effects, 
such as surface roughness scattering). However when the 
boundary is deformed (11) the angle of incidence of a ray 
associated with a lasing mode fluctuates and (at these 
deformations) is most likely to collide with the boundary 
below the critical angle of incidence at or near a loca- 
tion of high curvature (0 = 0°, 180°). Fig. 3B shows the 
calculated intensity pattern (the modulus squared of the 
electric field) for a typical "whispering-gallery" mode in 
a deformed cylinder laser with e = 0.06 (the calculational 
technique will be discussed below). The pattern shows 
clearly the enhanced emission intensity in the near-field 
in the vicinity of (<f> = 0°, 180°). The experiments are 
sensitive to the far-field intensity distribution, which de- 
pends also on the angle of refraction at the points of 
high curvature. Both the ray and wave calculations dis- 
cussed below indicate that at this deformation all "whis- 
pering gallery" modes with high output coupling have 
a minimum in emission intensity in the far-field around 
<t> = 0° and enhanced emission between 45° and 90°. The 
observed experimental intensity pattern has this general 
trend (see Fig. 3A) , but a fully angle- and mode- resolved 
measurement of the far-field pattern and a detailed com- 
parison with theory is difficult due to the generally low 
optical power and the many modes which contribute to 
the laser signal in this regime of deformations. (A de- 
tailed discussion of the spectral properties is given in the 
next section.) 

At higher deformations (e > 0.14) we detect a much 
stronger and qualitatively different directionality. Fig- 
ure 3C shows the actual angle-resolved far-field pattern 
(one quadrant) of one circular and two deformed lasers. 
For the laser displayed in Fig. 1, we obtain a power in- 
crease by a factor of 30 into an emission angle of <j> = 42° 
compared to <f> = 0°. The angular width of this direc- 
tional emission is w 23°. Around 0° we observe a clear 
minimum of the emission, and a smooth sloping plateau 
towards 90°. 

At these large deformations, a typical ray characteriz- 
ing a "whispering-gallery" mode escapes in less than 10 
collisions with the boundary as discussed in the theory 
section. This ray escape is now approximately isotropic 
and would seem unlikely to lead to the increased emission 
anisotropy observed experimentally. Because the general 
ray motion is furthermore highly chaotic in most of the 
phase-space, the only plausible scenario for generating 
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Figure 3: (A) Peak output power of different lasers as a function of deformation. The power is collected around 0° (open symbols) and 90° 
(filled symbols) with a width of the fixed aperture of 15°. Two independent sets of lasers are presented for each orientation of the aperture. 
Both curves rise approximately exponentially, as indicated by the dashed line-fit (26). Inset: Spectrum in linear scale obtained near peak 
optical power from a cylinder laser with low deformation (e ?s 0.04). The close mode spacing observed in the spectrum is a result of several 
lasing whispering-gallery-type modes. The displayed linewidth is limited by the experimental set-up and data acquisition system. (B) 
False color representation of the radiation intensity pattern of a chaotic "whispering-gallery" mode for a deformed cylinder with e = 0.06 
and length of the minor axis of 50 /im. Red indicates high intensity, dark blue minimum intensity on a logarithmic intensity scale. The 
computational technique is explained in the theory section of the main text. (C) Symbols indicate the measured angle-resolved far-field 
pattern (one quadrant) of a circular (A) and two deformed lasers with t = 0.14 (O) and e = 0.16 (•). The measurements presented 
here have been taken at a constant current level, at which the (deformed) lasers displayed pure single mode emission. However, the 
far-field shows qualitatively the same characteristic directionality at a current level corresponding to peak optical power. The data-sets are 
normalized to the value measured at zero degree. The data points are connected by splines (solid lines) for clarity. The dashed line is the 
far-field intensity pattern associated with the bow-tie mode shown in Fig. 3D, averaged over the experimental aperture. The calculation 
has been scaled to match the peak emission at 45° . The exact angular position of this maximum should be quite sensitive to the precise 
shape of the boundary near the "bow-tie" impact points, and we attribute the off-set between the measured and calculated peak positions 
primarily to the small deviation between our model and the actual shape and some uncertainty in the precise measurement of the angle. 
Furthermore, at present, we do not have a full understanding of the discrepancy in the intensities of the secondary peak between calculation 
and experiment. (Left inset) Logarithmic plot of the measured spectrum at maximum power (power P versus wavelength A) of a laser 
with e 0.16. Six equally spaced modes, with mode spacing AA = 40.4 nm, are observed. This mode separation is in good agreement 
with the value of 39.5 nm calculated for a "bow-tie" orbit corresponding to the calculated intensity pattern of Fig. 3D. (Right inset) 
The polar coordinate system is oriented such that = 0° indicates the direction along the elongated (major) axis, and <f> = 90° denotes 
the direction of the compressed (minor) axis. (D) False color representation of the intensity pattern of a "bow-tie" -mode for e = 0.15 and 
length of the minor axis of 50 fim. Note that the cross-over to the asymptotic far-field pattern of Fig. 3C (dashed line) is rather slow, and 
that certain features such as the modulated intensity at <j> = 90° vanish in the far-field. The color scale (red - high intensity, blue - low 
intensity) is unrelated to the color scale of Fig. 3B. 
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directional emission is for the lasing modes to be associ- 
ated with the small regions of stable, regular motion that 
still remain. 

For the range of deformations e « 0.12 to 0.23, there 
exist only two such regions. The first is the vicinity of 
the basic diametral orbit running along the minor axis 
of the resonator. The associated modes are the trans- 
verse modes of the stable, curved mirror Fabry-Perot 
resonator. However, these modes correspond to normal 
incidence at the boundary and, as such, would result in a 
peak emission at 90° in the far-field, much in contrast to 
the observation. Furthermore, the low reflectivity of the 
boundary at normal incidence combined with the short 
length of the minor axis result in a too high threshold for 
laser action. 

The second region is in the vicinity of the stable four- 
bounce periodic orbit with the shape of a "bow-tie" in 
real space. The intensity pattern of a representative 
"bow-tie" mode is shown in Fig. 3D. This orbit comes 
into existence by bifurcation from the diametral orbit at 
e « 0.10 and has four equal (in absolute value) angles 
of incidence on the boundary. At e w 0.12 this angle 
w 12.5° and it is well below the critical angle, but as the 
deformation increases to e rts 0.15, this angle increases to 
approximately the critical angle, \c ~ 17.5°. This change 
results in a sufficiently high reflectivity of the boundary 
to allow for laser action. For e of 0.125, 0.14, and 0.15, 
values of the reflectivity of 0.45, 0.59, and 0.72 are cal- 
culated, respectively. In fact, this increase in reflectivity 
with deformation should lead to a reduction of the laser 
threshold. 

When the radiation intensity pattern of a "bow-tie" 
mode is averaged according to the experimental condi- 
tions, we find reasonable agreement between the experi- 
mental and theoretical far-field directionality (Fig. 3C). 
We conclude that the laser emission at high deformations 
originates from the newly observed "bow-tie" modes. 
The spectral properties of the emission provide further 
confirmation of this fact, as discussed in the following 
section. 

The "bow-tie" orbit is just one of several orbits that 
move around the minor axis in a librational motion (that 
is without a fixed sense of rotation) as opposed to the 
rotational motion of conventional "whispering gallery" 
orbits. With higher index of refraction or different shape 
deformations, modes associated with other librational or- 
bits may be relevant to lasing; hence we will refer to the 
"bow-tie" as one of a class of "librational" modes. 

In general one would expect that the threshold current 
density Jth should have a minimum for e = (circular 
case), increase with deformation until e w 0.1, due to the 
increase in outcoupling loss, and then decrease due to the 
gradually increasing reflectivity of the "bow-tie" modes. 
In the range of e w 0.12 to 0.2, the observed decrease in 
Jth from « 5 kAcm~ 2 to 4 kAcm~ 2 is consistent with 
this expectation. However, in the "whispering gallery" 



range of deformations e = to e s=s 0.08 the measured 
decrease in J t h from « 7 kAcm~ 2 to 6 kAcm -2 is in 
contrast with the expectations. There are several issues 
that complicate the interpretation of the threshold data. 
First, there is a finite lateral current spreading resistance, 
effectively reducing the current density towards the edge 
of the disk outside the contact region. Second, the mode 
confinement factor within the active region is expected 
to be reduced in the outermost parts of the waveguide 
due to the true three-dimensional nature of the waveg- 
uide, increasing locally the threshold current density. As 
a consequence, the actual threshold current density for 
a given mode depends on its spatial distribution within 
the resonator. Our experimental Jth, on the contrary, 
is always calculated by dividing the value of the thresh- 
old current by the geometrical area of the actual device. 
Finally lowering of the threshold current density will in 
general also lead to an increase in the maximum output 
power due to the higher available range of drive currents. 

The threshold current density of the QC-laser is given 
by Jth = {oi wav + a out )/gr, with a wav being the waveg- 
uide loss, gr the average modal gain coefficient, and 
otout, the outcoupling loss, which strongly depends on 
the distributed reflectivity of the boundary, which in 
turn depends strongly on the deformation e, and effec- 
tive length of the resonator. From the laser threshold 
and the computed value of gr (6.72 at 10 _3 cmA~ ), 
the quality-factor ("Q- value") can be calculated as Q — 
{2tt n)/[\ (a wav + aw)], where n is the effective refrac- 
tive index and A the wavelength. The above cited thresh- 
old current densities then result in Q-values ranging from 
w 850 to 1500. It should be noted that waveguide losses 
are usually dominant in QC lasers due to the high doping 
levels which increase free carrier absorption. 

Spectral properties. Besides the increase in output 
power and directionality, the increasing deformation also 
influences the spectral properties of the lasers. These 
further confirm the existence of two different regimes, 
manifested so clearly in the different types of far-field 
patterns. 

At low deformations we obtain a complex, dense modal 
spectrum. The lasers are multiple-mode starting from 
threshold, with close mode spacings, and show up to 
10 almost equally strong modes at the maximum opti- 
cal power (see inset to Fig. 3A). This close mode spacing 
cannot be understood from one fundamental set of longi- 
tudinal (azimuthal) whispering gallery modes only since 
fitting an integer number of wavelengths along a single 
closed ray path would result in a regular comb of modes 
with significantly larger spacings. We therefore attribute 
the spectrum to the lasing of several different longitudi- 
nal ( "azimuthal" ) and transverse ( "radial" ) modes. 

At large deformations the lasers are single-mode until 
approximately twice the threshold current, and show at 
most two to three strong modes at maximum power. The 
onset of additional modes is accompanied by a kink in 
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the light output-current characteristic; one can be seen 
in Fig. 2B. The cross-over between the two spectrally 
characteristic regimes again occurs around e w 0.12. 

The multiple mode behavior of the highly deformed 
lasers is consistent with the emission from "bow-tie" 
modes. The logarithmic plot of a spectrum in this regime 
(inset Fig. 3C) reveals six equally spaced modes, with 
mode spacing AA = 40.4 nm. The expected theoretical 
value is calculated assuming that adjacent modes differ 
by one wavelength along the path-length of the "bow- 
tie". This analysis yields a spacing of 39.5 nm, in ex- 
cellent agreement with the experiment, considering the 
uncertainty in the effective refractive index. 

The "bow-tie" modes can easily be distinguished from 
transverse modes of the diametral curved mirror Fabry- 
Perot resonator along the minor laser axis (length L). As 
noted above, they originate from a period-doubling bifur- 
cation of the latter, as will be discussed below in more 
detail, leading to approximately twice the optical path 
length. As such, the "bow-tie" mode spectrum displays 
approximately half the mode spacing one would expect of 
the standard Fabry-Perot modes, AA = X 2 /(2nL) w 82 
nm. 

To summarize, the experimental data show that impos- 
ing a flattened quadrupolar deformation onto semicon- 
ductor micro-lasers substantially improves their power 
output and directionality. In the favorable directions of 
the far-field a power increase of up to three orders of 
magnitude was obtained. This dramatic result could be 
achieved by exploiting the complex ray dynamics - first 
for chaotic "whispering-gallery" modes, then for "bow- 
tie" modes - of the deformed resonators. An in-depth 
theoretical discussion of the subject is given in the fol- 
lowing section. 

Theory. The intensity patterns shown in Fig. 3, 
B and D, were obtained by numerical solution of the 
Helmholtz equation for the TM polarization resonances 
at A w 5.2 ^m of a deformed dielectric cylinder with the 
dimensions and index of refraction (n = 3.3) correspond- 
ing to the experimental structures. These solutions are 
obtained by matching the internal and external electric 
fields and their derivatives at the surface of the semicon- 
ductor, along with the additional constraint that there 
is no incoming wave from infinity. The latter constraint 
implies that the wavevector must be complex, with the 
imaginary part giving the decay rate or Q-value of the 
resonance (28). 

To obtain a full theoretical understanding of these res- 
onances it is helpful to divide the problem in two parts. 
First, we consider the properties of the "bound-states" of 
the system, corresponding to the discrete solutions which 
would exist if the cavity were completely closed and the 
electric field were zero outside the cavity. Then we must 
understand how these states are altered by the possibility 
of escape to infinity by refraction. 

The first point is precisely the issue of understanding 



the solutions of the wave equation within a "billiard" . 
This problem corresponds to a resonator with a mirror 
reflectivity exactly equal to unity. With these "hard- 
wall" boundary conditions, the Helmholtz wave equation 
is identical to the Schrodingcr equation of quantum me- 
chanics. 

When the cross-section of the cylinder is deformed 
from circularity, the wave equation is no longer separa- 
ble into three ID differential equations and the solutions 
in the plane transverse to the cylinder axis are no longer 
specified by pairs of quantum numbers (or mode indices). 

One can still obtain a numerical solution, that is by 
representing the solution in a large basis set of states 
and diagonalizing the resulting matrix equations. How- 
ever, if this approach is used alone it is difficult to extract 
any physical understanding of the bound states or "true" 
resonances, now taking into account the electric field out- 
side the resonator. In fact, the solutions shown in Figs. 3, 
B and D were predicted first by a completely different 
theoretical approach, before they were found by numer- 
ical search. This different approach, which has been pi- 
oneered in physics (29-31) and physical chemistry (32) 
during the past two decades, is to study the short wave- 
length limit of the problem (ray optics for the Helmholtz 
equation, Newtonian mechanics for Schrodinger's equa- 
tion) and try to develop a systematic understanding with 
semiclassical methods. The use of semiclassical methods 
is justified in our system because the wavelength of light 
in the material (« 1.6 /xm) is much smaller than any of 
the geometric features of the resonators. Moreover, stan- 
dard perturbation techniques are not applicable because 
the deformation causes a shift in the resonance frequen- 
cies that is large compared to the resonance spacing. 

When the optical wave equation is non-separable, the 
corresponding ray motion typically exhibits fully or par- 
tially chaotic dynamics, just as the classical limit of 
a non-separable Schrodinger equation typically gives a 
chaotic classical mechanics; this sub-field has become 
known as "quantum or wave chaos theory." The sta- 
tionary states of these so-called "quantum billiards" have 
been studied extensively in this context. Here we will dis- 
cuss the ray-optics properties of the billiards correspond- 
ing to the laser resonators studied above, with the goal 
of understanding the cross-over between emission from 
"whispering-gallery" to "bow-tie" modes which occurs in 
this system. 

The relevant billiards are smooth deformations of the 
circular billiard. Initially we neglect the possibility of 
escape. Rays will simply propagate indefinitely within 
the billiard, satisfying the law of specular reflection at 
collisions with the boundary. When the circle is unde- 
formed, angular momentum is conserved in this motion. 
The angle of incidence, Xi is the same at each collision, 
and the orbit traces out an annulus bounded by a circu- 
lar caustic of radius R sin(x), where R is the radius of 
the circle. The corresponding wave solutions are the ordi- 
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nary Bessel-functions indexed by the angular momentum 
quantum number. 

When the boundary conditions are changed to include 
refraction, then rays incident with \ greater than the 
critical value Xc given by sin(xc) = l/ n wm remain 
trapped by total internal reflection, whereas rays with 
sin(x) < l/n will rapidly escape by refraction according 
to Snell's law. 

To illustrate the circular and deformed case in a unified 
manner we represent the ray motion in phase space using 
the surface of section (SOS) method (31,33), in which 
every time a ray collides with the boundary both the 
azimuthal angle (0), at which it hits, and its angle of 
incidence (x) with respect to the boundary are recorded. 
Following an ensemble of hundred trajectories for 200 
bounces then gives a good picture of the global dynamics 
in phase space. 

The generic behavior of smoothly deformed circular 
billiards in this representation is shown in Fig. 4, where 
again we neglect the possibility of escape in calculating 
the SOS. For the circle (Fig. 4A) the SOS is trivial, each 
trajectory gives a straight line corresponding to the con- 
served value of sin(x), except for trajectories with a chord 
angle (2%) equal to a rational fraction, p/q, of 2tt. Such 
trajectories will close after q bounces and are referred 
to as "period-g" orbits. All such orbits in the circle 
are marginally stable and exist in infinite families cor- 
responding to arbitrary rotations of any one orbit in the 
family. Several period-2, period-3, and pcriod-4 orbits 
are indicated in the SOS of Fig. 4A; the period-2 orbits, 
which are very important in the discussion below, just 
traverse the diameter of the circle. 

In all of the SOSs in Fig. 4 we have indicated in red 
the horizontal line corresponding to the critical angle, 
sin(x c ) = l/n = 0.30. Trajectories which fall below that 
line in the closed billiard will escape from the semicon- 
ductor. Trajectories above stay "forever" trapped within 
the resonator [in this approximation, which neglects weak 
evanescent leakage or tunneling of photons (34)]. When 
the circle is deformed the ray dynamics in the billiard un- 
dergoes a transition to partially chaotic motion. If the de- 
formation is smooth and the curvature of the boundary is 
always convex, it can be shown rigorously that the phase- 
space still has non-chaotic "whispering-gallery" modes 
for values of sin(x) sufficiently close to one (35). 

The specific form of the deformation is unimportant for 
the qualitative physics; we use the flattened quadrupolar 
deformation, which describes well the experiment. One 
sees the effect of a deformation of e = 0.06 in Fig. 4B. For 
sin(x) > 0.7, there remain many unbroken (continuous) 
curves traversing the full surface of section which corre- 
spond to "whispering-gallery" modes that survive only 
slightly deformed from the circle. These are "whispering 
gallery" orbits of the familiar type, which are confined 
near the rim of the resonator, have a true caustic and will 
circulate in one sense indefinitely. However one also now 



sees the signature of isolated stable and unstable periodic 
orbits in the motion. The deformation destroys the infi- 
nite number of periodic orbits in each family and leaves 
just an equal number of stable and unstable orbits. The 
stable orbits are surrounded by closed curves ( "islands" ) 
which indicate the oscillatory motion of nearby trajecto- 
ries around the stable periodic orbit. The simplest ex- 
ample in Fig. 4B are the two islands around the stable 
(short) diametral orbit which collides with sin(x) = 
at <p = 90° • The unstable orbits generate regions of 
chaotic motion nearby them, which correspond to the 
grainy structureless regions of the SOS. The most visible 
example in Fig. 4B extends around the period-2 islands, 
reaching the sin(x) = axis at the location of the un- 
stable (long) diametral orbit (which has sin(x) = and 
= 0, 180°). The "bow-tie" modes which we have fo- 
cused on in the previous sections would correspond to a 
four-bounce orbit centered on the diametral orbit around 
<f) = 90°, but no such orbit exists at this low deformation. 

To confirm that the relevant resonances at this low 
deformation arc of the "whispering gallery" type, one 
can generate a phase-space representation of the inten- 
sity pattern of Fig. 3B, called the Husimi function(35). 
For the resonance with deformation e = 0.06 shown in 
Fig. 3B, this function (shown in Fig. 4E) demonstrates 
that the ray motion corresponding to this state is spread 
out in the large chaotic region just mentioned. Note that 
since the chaotic region extends through sin(x) = an 
orbit in this region of phase space will eventually change 
its sense of rotation and is not a "whispering gallery" 
orbit in the familiar sense. However the Husimi func- 
tion of Fig. 4E does not have support near sin(x) = 
indicating that escape occurs before this reversal of cir- 
culation can happen; hence the corresponding real-space 
intensity pattern (Fig. 3B) does have a minimum in the 
center bounded by an approximate caustic. Note fur- 
ther that this orbit lies entirely outside the influence of 
the central diametral orbit and collides with all regions 
at the boundary, thus it may reasonably be termed a 
"chaotic whispering gallery" orbit. 

At a deformation e = 0.10, the "bow-tie" orbit ap- 
pears at a period-doubling bifurcation (33) of the stable 
diametral orbit. In this case, it is a non-generic period- 
doubling bifurcation (36) in which a new stable orbit of 
twice the period is born (the "bow-tie"), while simul- 
taneously two new unstable, V-shaped, period-2 orbits 
("birds") are born (see Fig. 5). Such period-doubling 
bifurcations are well-understood and can be described 
quantitatively within the general formalism of non-linear 
hamiltonian dynamics (33). However, in this case one 
can also use a more elementary argument familiar from 
resonator theory. The stable (vertical) diametral orbit 
supports standard Gaussian Fabry-Perot modes which 
are too low-Q to lase in this structure, due to the rel- 
atively low reflectivity at normal incidence. When the 
radii of curvature at the two contact points of this orbit 
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Figure 4: (A to D) Poincare surface of section representing the motion of an ensemble of rays in phase space for the flattened quadrupolar 
billiard, neglecting the possibility of refractive escape. Regions of stable or regular motion are indicated in green, and regions of chaotic 
motion are indicated in blue. is the polar angle in radian as defined in Fig. 3C, right inset. (A) The undeformed (circular) cylinder. Each 
trajectory collides with the boundary at a fixed value of the angle of incidence, sin(x), sometimes closing on itself and forming a periodic 
orbit, otherwise passing arbitrarily close to any point on the boundary and forming a line in the surface of section. Several members of the 
infinite families of period-2 (circles), period-3 (triangles), and period-4 (squares, diamonds) orbits are shown. These members survive to 
nonzero deformation, filled symbols represent orbits that will give stable islands, empty symbols those that will be unstable and generate 
regions of chaos. The red line represents the escape condition, sin(x) = 1/n; in the true resonator rays below that line would rapidly 
escape by Fresnel scattering. (B) The phase space for deformation e = 0.06, corresponding to the calculation of Fig. 3B. The two major 
islands at polar angle <j> = ±7r/2 (±90°) correspond to motion around the stable diametral orbit. Just above these islands is the chaotic 
region generated by the unstable diametral orbit at cf> = ±7r(±180°). This region contains the chaotic "whispering gallery" modes. (C) 
For e = 0.125, somewhat after the bifurcation of the diametral orbit (at e = 0.10) that gives rise to the "bow-tie" orbits which are seen 
clearly as the four islands at sin(x) = 0.22. These islands are sufficiently below the critical line that the corresponding modes would be too 
short-lived to lase. Note that there are four symmetric islands for negative sin(x) which are not shown; any one "bow-tie" orbit only visits 
four of the islands, two with positive and two with negative sin(x), but the same path is traced in either case. (D) The behavior at e = 0.15 
corresponding to the data and calculation of Fig. 3, C and D. Now the "bow-tie" islands have moved up to the critical line, increasing the 
lifetime of the corresponding modes and allowing them to lase. (E) Husimi function corresponding to the resonance of Fig. 3B at e = 0.06; 
this function clearly represents a chaotic "whispering gallery" state localized in the chaotic region. The Husimi function translates the 
real-space electric field intensity pattern into a probability density in phase-space. The resulting function is illustrated by a color scale, 
where red is high intensity. (F) The relation between the highly directional resonator mode shown in Fig. 3D and the islands shown 
in Fig. 4D can be demonstrated by means of the Husimi function shown in Fig. 4F. This function is centered on the "bow-tie" islands. 
The minimum at the very center indicates that this "bow-tie" mode has an oscillatory motion transverse to the "bow-tie" path; this is 
consistent with the intensity pattern of Fig. 3D which exhibits four transverse oscillations. 
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Figure 5: The three orbits born at the period-doubling bifurca- 
tion of the stable diametral orbit, the stable "bow-tie" (red) and 
the two unstable V-shaped "birds" (blue, green). This occurs at 
a deformation e = 0.10. The birds, being unstable, will not gen- 
erate long-lived resonator modes; however, the stable "bow-tie" 
generates modes with directional properties and spectral spacing 
in excellent agreement with the experiment as discussed above. A 
key feature of the "bow-tie" is that it does not exist until the res- 
onator is substantially deformed so that the confocal condition is 
reached for the stable diametral orbit, as discussed in the text. 

become equal to the distance between them (the minor 
axis L) we reach the confocal condition (37) at which 
marginally stable families of "bow-tie" and V-shaped or- 
bits, all of length 4L, come into existence. For the flat- 
tened quadrupole this occurs at e = 0.10. For slightly 
larger deformations these orbits leave the vicinity of the 
diametral orbit and do not correspond to small deforma- 
tions of diametral orbits. Such orbits are not typically 
discussed in Fabry-Perot theory (38). But here, because 
the boundary creates a full 180 degree "mirror" with a re- 
flectivity which increases at oblique incidence, the modes 
associated with the remaining stable "bow-tie" orbit are 
higher-Q than the simple Fabry-Perot modes and can lase 
when the latter do not. Because they require a doubling 
of the radius of curvature at the minor axis, they do not 
exist at small deformations. 

The SOS for e = 0.125 shown in Fig. 4C is taken just 
after the bifurcation of the diametral orbit, showing the 
emergence of the stable "bow-tie" , which has the feature 
that its angle of incidence is the same for all four bounces 
(see caption to Fig. 4C). Note, however, that for this 
deformation sin(x) = 0.22, which is still well below the 
critical angle. 

When the deformation is further increased to e > 0.14 
the "bow-tie" orbit has moved upward in the SOS so that 
it is centered near the critical angle (see Fig. 4D). The 
reflectivity of the corresponding modes will increase to 
a value comparable to the "whispering-gallery" modes 
at the same sin(x) in the circle. Therefore we expect 
a turn-on of the laser emission from this mode. Note 
that the "bow-tie" orbit now represents the only large 
stable island at or above the critical angle in the SOS. 
Moreover, it is now well separated in phase space by a 
chaotic region from the fundamental diametral orbit from 
which it originated. The plot in Fig. 3D shows the high- 
intensity regions concentrated on this orbit. In Fig. 4F, 
we show the phase-space projection of this mode, which is 
concentrated in the vicinity of the islands corresponding 



to the "bow-tie" orbit. 

The modes corresponding to the "bow-tie" orbit are 
not simply higher order transverse Fabry-Perot modes; 
the latter would correspond to quantized oscillations 
within the island around the diametral orbit. Moreover, 
as noted above, the "bow-tie" orbit is rather different 
from the "whispering-gallery" orbits because the sense 
of rotation of the "bow-tie" orbit is not constant; it rep- 
resents a librational rather than a rotational motion. 

The existence and stability of the "bow-tie" orbit is 
relatively insensitive to the precise shape of the bound- 
ary, so we expect these modes to be generic to deformed 
cylindrical resonators. For the flattened quadrupole, the 
stable "bow-tie" exists in the range of deformations from 
e = 0.10 to e « 0.23. Its directions of peak emission, 
though, are sensitive to the precise shape of the res- 
onator; the degree of sensitivity will be the subject of 
future studies (38). Nevertheless, reasonable agreement 
between theory and experiment has been obtained for 
the far-field directionality using the flattened quadrupo- 
lar shape (Fig. 3C). 

As noted above, for the range of deformations at which 
the stable "bow-tie" orbit exists, it represents the only 
substantial islands of stability in the region of phase- 
space close to the critical value of total internal reflec- 
tion; thus it is difficult to find a competitive mechanism 
for the highly directional modes we observe. At lower 
deformations, other librational modes exist and may be 
important in the cross-over from "whispering gallery" to 
"bow-tie" emission. 

Highly directional emission from low refractive index 
resonators was discussed in earlier theoretical work by 
several of the authors, and tested in experiments on las- 
ing dye-jets (12). However, the origin of directionality 
at high deformations in the high refractive index res- 
onators discussed in the present paper is qualitatively dif- 
ferent from the mechanism studied in this earlier work. 
In resonators with indices of refraction n < 2, the es- 
cape line corresponding to sin(x) = 1/n is much higher 
in the surface of section. Therefore a ray escaping 
from a "whispering-gallery" mode must traverse a much 
smaller fraction of the chaotic sea to escape. It has been 
shown (11,12) that in this case the motion is not pseudo- 
random, and highly directional emission from near the 
points of highest curvature results. However, in the high- 
index materials of the present work, it is necessary to 
reach much lower angles of incidence within the resonator 
to escape and we now find that the escape direction for 
rays starting far from the critical angle is effectively ran- 
dom, at least for the deformations where the "bow-tie" 
orbit exists. This is demonstrated by the chaotic scatter- 
ing map shown in Fig. 6. As is explained in the caption 
to Fig. 6, this map suggests strongly that highly direc- 
tional modes of the "whispering-gallery" type are not 
easily achieved at high deformations in such resonators 
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Figure 6: Color representation of the directionality of escaping rays 
for the phase space of the flattened quadrupole at e = 0.15. Blue 
represents initial conditions leading to escape into the far-field at 
0° and red at 90° . Note the clear demarcation between a pseudo- 
random region with rapidly fluctuating escape direction for initial 
angles sin(x) > 0.70 and a regular region where the escape di- 
rection varies smoothly and relatively slowly. In the regular re- 
gion, which primarily corresponds to librational motion, escape is 
so rapid that chaos cannot fully develop; in contrast, initial condi- 
tions in the "whispering gallery" region above sin(x) = 0.70 must 
traverse more of the chaotic sea and cannot generate highly direc- 
tional escape. This behavior is completely different from low-index 
resonators with the same magnitude of deformation (11). The reg- 
ular region can generate directional emission, but only for states 
localized by islands, such as the "bow-tie" states. 

made from semiconductor materials, although such 
modes exist and dominate the lasing properties at the 
same range of deformations for lower index materials. 
Conversely, modes such as the "bow-tie" resonance, 
which are related to librational orbits, all reside well be- 
low sin(x) = 0.5 and as such would experience too little 
reflectivity from the boundary to reach laser threshold 
in low-index materials. Note finally that the "bow-tie" 
modes are confined away from the points of highest curva- 
ture in the resonator and thus display a minimum in the 
near-field intensity at these points (Fig. 3D) in contrast 
to the "whispering gallery" modes (Fig. 3B) which have 
high intensity in the near-field at these points. There- 
fore the two types of modes should be easily distinguish- 
able if near-field measurements could be made. This is 
a demanding task for lasers in the mid-infrared region of 
the spectrum and will be the subject of future research. 
Finally, it should be emphasized that there is no funda- 
mental reason that such resonators should not be equally 
effective as microcavities at visible wavelengths. 
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